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ABSTRACT 


This paper will detail the navigation and dispersion design and analysis of the first Orion crewed 
mission. The optical navigation measurement model will be described. The vehicle noise includes 
the residual acceleration from attitude deadbanding, attitude maneuvers, C'O 2 venting, waste- 
water venting, ammonia sublimator venting and solar radiation pressure. The maneuver execution 
errors account for the contribution of accelerometer scale-factor on the accuracy of the maneuver 
execution. Linear covariance techniques are used to obtain the navigation errors and the trajectory 
dispersions as well as the DV performance. Particular attention will be paid to the accuracy of the 


delivery at Earth Entry Interface and at the Lunar Flyby. 


The Orion Exploration Mission-2 (EM-2) will be the first Orion mission to carry crew and there- 
fore the mission will be subject to high level of scrutiny. The Orion vehicle is thereby required to 
return the crew safely in the case of loss of communication with the ground. As such, it needs to 
be able to navigate autonomously, independent of ground-based measurements, utilizing on-board 
sensors subject to stringent mass/power/volume constraints. Since the vehicle will be carrying opti- 
cal cameras, the cislunar navigation system is designed to use images obtained from these cameras, 


in particular planetary limb measurements and therefore presents unique presents challenges. 
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This paper details the navigation design and analysis of the Orion Exploration Mission 2 (EM-2). 
The characteristics of a Multi-Trans-Lunar-Injection trajectory with an accompanying free-return 
trajectory present unique GNC challenges that will be addressed in the paper. This is depicted in 
Figure 1. Whereas the Orion sensor complement includes two star trackers, the chosen star trackers 
have a very limited field-of-view. As such, they don’t lend themselves to cislunar optical naviga- 
tion, which needs fields-of-view in excess of 20 degrees. A dedicated optical navigation camera 
is mounted on the same optical platform as the star trackers. The paper will present analysis (and 
resulting design) that has been performed to understand the performance of the Orion system for 


EM-?2, with particular attention paid to entry flight path angle constraints and the DV performance. 
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Figure 1. The Orion EM-2 Trajectory Profile 


The vehicle noise model includes five sources: attitude deadband residual accelerations, attitude 
maneuver residual accelerations, pressure swing adsorption (PSA) vents which are used to vent 
the CO, in the crew module, waste water vent accelerations and ammonia sublimator vents. For 
EM-?2 and following, the Pressure Swing Adsorption (PSA) cycles automatically approximately 
every 10 minutes, less during crew sleep periods, more when the crew is awake — the analysis 
will also include periods of astronaut exercise each day. The waste water vents are expected to 
occur every 8 hours. The ammonia sublimator vents are expected to only occur 20 minutes prior to 


entry. This analysis will assume a 24-jet thruster configuration. Each of the these errors has been 
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modeled from first-principles and since it is not known a priori when these events will occur, they 


are spread out over the entire mission. The paper will discuss this in detail. 


The maneuvers are modeled in terms of scale-factor, bias, misalignment, and noise. Long ma- 
neuvers present a particular challenge in meeting accuracy requirements. The accelerometer bias, 
in particular, is used to model the maneuver scale factor errors which are the major component of 


maneuver execution error. The paper will discuss this in detail. 


The paper will detail the optical camera measurement model used to perform the navigation. 
The measurements are the planetary centroid and the distance to the planet. This is obtained 
from the onboard image processing software. The onboard navigation filter then processes these 


measurements to estimate the inertial position and velocity of the Orion vehicle. 


In an emergency situation, during a loss of communication scenario, the primary objective is to 
secure the safety of the crew. This subsequently translates into a flight-path angle requirement at 
entry interface (EI) for a direct entry. A direct entry, as opposed to a skip entry, reduces the risk of 
the capsule bouncing back into space, and allows for a greater margin on the flight-path angle at 


EI. 


Linear covariance analysis will be used to perform the navigation and dispersion analysis for the 
EM-2 mission. Particular attention will be placed to the trajectory dispersions at the lunar flyby 


and at entry interface. 


The accuracy of the flight-path angle at EI is driven by several factors including the navigation, 
targeting, and burn execution errors at the time of the last mid-course maneuver, and unaccounted 
trajectory perturbations between the last mid-course maneuver and EI. Apollo missions tolerated 
a maximum flight path angle error at EI of +1 degree, with half of this error allocated to naviga- 
tion. A similar criterion is employed in this study. However, for Orion, the EI dispersion error 
requirement with optical navigation is 0.3 degrees (30). The paper will evaluate the accuracy of 


the delivery in the presence of the aforementioned errors. 
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LINEAR COVARIANCE ANALYSIS 


This investigation is performed using linear covariance (LinCov) analysis techniques.''? The 
state vector is 


.a4n- ve ot m,, mem ete (1) 


where r is the inertial position of the vehicle with respect to the primary body, v is the correspond- 
ing inertial velocity with respect to the primary, @ is the attitude error, m,, is the misalignment 
of the optical instrument, b;, is the bias position of the instrument with respect to the navigation 
base, beenz is the bias of the planetary centroid measurement, and b,, is the bias of the planetary 


diameter measurement. 


Neither the attitude error 6@ nor its uncertainty are integrated in this analysis. The nominal 
attitude is known at any time and it does not need to be calculated. The attitude estimation error 
covariance 1s treated as a constant and is driven by the star tracker accuracy. The attitude navigation 
dispersion covariance is also treated as constant and is given by the attitude control dead-band. Be- 
fore the star elevation is determined, the vehicle slews in preparation for measurement acquisition. 
This attitude maneuver is performed by the onboard thrusters and is assumed to be instantaneous. 
Due to thruster misalignment, this maneuver adds uncertainty to the translational states. After the 
batch of measurements is available, the vehicle returns to its nominal attitude. In linear covariance 
analysis, the difference between the true state and the nominal state is defined as the environment 
dispersion 


6x =x—X. (2) 


The difference between the estimated state and the nominal state is defined as the navigation dis- 


persion 


Ox £X—. (3) 


Finally, the difference between the true state and the estimated state, is defined as the estimation 
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error, sometimes referred to as the onboard navigation error 


e=x-x. (4) 


Following the standard Kalman filter assumptions, the difference between the nominal and esti- 


mated models is represented with zero-mean, white noise. The estimated state evolves as 


x = f(x), (5) 


where f is a nonlinear function representing the system dynamics as modeled by the filter. The 


evolution of the nominal state is modeled as 


x = F(X) =£(%) +0, (6) 


where f is a nonlinear function representing the state dynamics as modeled in designing the nom- 
inal trajectory. The nominal dynamics f may be higher fidelity than the filter’s dynamics f. The 
vector vu represents the dynamics modeled in the nominal trajectory but neglected in the filter mod- 
els. In Kalman filtering, the difference between the true dynamics and the filter’s dynamics is called 
process noise. While these unmodeled dynamics are not actually white noise, they are modeled as 
such. The power spectral density of process noise is then tuned to achieve good performance. The 
same procedure is used here. In order to capture the difference between the two dynamical models, 
v is modeled as a zero-mean white noise process with power spectral density Q. The goal is to 
represent the increased value of the navigation dispersion during propagation due to the difference 


between the nominal and filter’s dynamical models. 


The evolution of the navigation dispersion can be approximated to first-order as 


6x =X —X = £(X + 6k) — £(X) —v © F(X)6x— v. (7) 


The evolution of the navigation dispersion covariance is governed by 
P = F(x)P + PF(x)" +Q. (8) 
Similarly, the true state is modeled to evolve as 
x=f(x)+v. (9) 
The evolution of the estimation error is given by 
@=x—x~ f(x) + F(x)(x — x) +v — f(x) — F(x)(& —x) = F(X)e+v. (10) 


Vector v is modeled as zero mean white noise with power spectral density Q. The onboard covari- 
ance P evolves as 


P = F(x)P + PF(x) + Q. (11) 


Notice that the Jacobian F could be evaluated at the estimated state x instead of the nominal state 


X, as in the extended Kalman filter. 


Finally 
Ox =x—-xX 2 F(xX)bx+v—v (12) 


and P evolves as 


P = F(x)P + PF(x) + Q. (13) 


Notice that Q = Q + Q if v and v are assumed to be uncorrelated. 
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OPTICAL NAVIGATION MODEL 


Begin with the equations 


Yh tana, +b, +, = «/z+b,n+m (14) 


Yy = tanathtm = y/zthtm (15) 


p = VETPEF = [Prin —*5l (16) 


The first two, y;, and y,, are the centroid measurements and the last one is the range, which will be 
corrupted a bit further in this section by bias and noise to produce a range measurement. 


If the planet is perfectly spherical, the range to the planet would be given by 
(17) 


where RF, is the planetary radius and A is the angular diameter. The model assumed in this analysis 


is a pinhole camera model. The equivalent expression at the focal plane is 
tan— = — (18) 


where nq is the number of pixels that constitute the angular diameter and s is the pixel pitch (in 


units of pixels per length, assuming square pixels). This is seen in Figure 2. For what follows, we 


Figure 2. Angular Diameter Measurement 


define 


(19) 


where 0 is a non-dimensional quantity and is the “raw” measurement of the planetary diameter. 


Consistent with the centroid measurements. Given that the measurements are corrupted with noise, 


the actual planetary diameter measurement is 
S A 
6 = 6+b¢t+n = 2tan 5 + ba + na 


Given that Eq. (18) can be expressed as 


using elementary trigonometry, the above equation (Eq. (21)) can be written as 


6 
sin = 


2 V62 +4 
Eq. (17) and Eq. (22) are equated to solve for 6 yielding 


2R, 


or equivalently p can be written in terms of 6 as 


2 2 


and acknowledging the fact that the measurement is corrupted with noise and bias, we get 


9) 2 
= Reals 
Yo P Go) 


(20) 


(21) 


(22) 


(23) 


(24) 


(25) 


The Measurement Error Calculations 


In LinCov we model the pixel centroid measurements as 


y=] a — fsteno, =.= FSU (26) 
C= Ve f an =U, = fF Sys (27) 
so that 
du = —(fs)dyp (28) 
du = —(fs)dyp (29) 
With this in hand, 
Ou = (fs) Oy, (30) 
Oy = (fs) Oy, (31) 
and 
1 
Oy, = fs Ca (32) 
1 
Oy, = fs Oy (33) 


We have been assuming in LinCov that the pixel pitch is 4.8 x 10~° m/pixel and the focal length 


is 35.1 mm so that 


1 
a= 1.3675 x 1074 pixels! (34) 
s 


which is equivalent to 28.2 arc-seconds per pixel (= 1.3675 x 10~* - 3600 
DADS 10 
0, = 0.06+ so pixels 
p 
2 123¢ TO": : 
Oy = 0.06+ ae pixels 
p 


where p is the distance of the vehicle with respect to the planet (in ft). 


Therefore, the required accuracy on the two centroid ‘tangent angles’ is 


2.899 x 108 
Cy, = 8.205 x 10°°+ a 
: 108 
oy, = 8.205 x 10-6 - IO 
uv p /2 
Given that 
= fs 
the error in these are found to be 
1 
dé = —d 
fs 4 


and 


pa Rds ep —Redna 


d — 
C p “) p Na 


The standard deviation can be approximated as 


2 2 
— p = fy he 
Pp Nd 


Op d 


10 


- 180/7) and 


(35) 


(36) 


(37) 


(38) 


(39) 


(40) 


(41) 


(42) 


which for p >> R, becomes 


p 
=P. 43 
aes aa 9 (43) 
For the LinCov Analysis, we are using ¢,,, = 0.06 + 2.12 x 10'?/p pixel so 
0.06 O19 10 
fae (44) 


Na png 


where 74 is the expected planetary diameter in pixels. 


THE VEHICLE PROCESS NOISE 
Attitude Dead-banding Maneuver Errors 


Given a 24-jet ESA-SM config, with attitude dead-bands with 0.028 second on-time, a recent 


analysis recommended a process noise of 


O st Deadband — 34019 * 10-"° ft? /sec? (45) 


Attitude Slew Maneuver Errors 
For attitude maneuvers, over a period of 3.2 hours, a Monte Carlo analysis resulted in a process 
noise of 


OARS elias pe = 1 18 10 Ese" (46) 


To account for varying number of slews (= 2nNgaizevents) Over a period of time (T},ou;5), the process 


noise for attitude events with a 24-jet configuration is 


6.4 att events 
@ Aisiews = 1.5153 x 19710 (SSueets | ft? /sec3 (47) 


hours 
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Event Duration Between Vents | Cycles Per Day | Process Noise 
(minutes) (ft? /sec?) 

Sleep PSA 30 16.992 1390-10" 

Awake PSA ibs) 45.024 2.779 x 10-° 

Active PSA i Fes) 76.464 1251 slo. 


Table 1. Process Noise Allocations for Astronaut Activity 


PSA Puffs 


For EM-2 and following, the Pressure Swing Adsorption (PSA) cycles automatically approxi- 
mately every 10 minutes, less during crew sleep periods, more when the crew is awake. The PSA 
system can be commanded to force a complete desaturation (all vents open into space) and then 
they would likely not open automatically for 40-60 minutes. This would be needed during ground 
tracking passes and optical navigation passes. 

The detailed assumptions are as follows: the maximum impulse per vent is expected to be 4.52 
Ibf-sec. This can be further decomposed as 0.92 lbf-sec of initial impulse and 3.60 Ibf-sec of 
decay impulse (for a total impulse of 4.52 Ibf-sec). Assuming a mass of 1650 slugs, the AV per 
vent is 2.739 x 1073 ft/s. Assuming no preference for a particular direction, the per-axis AV is 
1.582% 107? ft/s. 

Assuming that each day the astronauts will be asleep for 8.5 hours (35.4 %), awake for 11.25 
hours (46.9 %) and in vigorous activity (i.e exercising) for 4.25 hours (17.7 %), the process noise 


allocations are expressed in Table 1. 


The Process Noise can be calculated as 


Qpsa = 35.4% * Qsteep + 46.9% * Qawake + 17.7% * Q Active = 4.010 x 107° ft?/sec® (48) 


The total number of vents per day is 138.48 which means that the ‘average’ frequency of the 


vents is 623.9 seconds or 10.4 minutes. 
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Waste Water Vents 


From EM-2 onward, crewed missions will also be required to perform waste water (urine) vents. 
Each of these vent events have a maximum impluse (per vent) of 10.2 Ibf-sec. This results ina AV 


per vent of 6.182 x 10~° ft/sec or 3.569 x 107° ft/sec per axis. The resulting process noise model 


is described in Table 2. 
Event Duration Between Vents | Cycles Per Day | Process Noise 
(hours) (ft? /sec?) 
Urine Vent 3 8 1.179 x 107° 


Table 2. Process Noise for Waste Water Vents 


Solar Radiation Pressure 


A simple model presented here, based upon Wertz,” is 


Fsrp = CrpsrA (49) 


where Cp is the coefficient of reflectivity (a value which varies between 1.0 (for perfect absorption) 
and 2.0 (for perfect reflectivity)), psp is a constant that contains the average/expected solar radia- 
tion pressure at 1 a.u. (in N/m?) and A is the expected exposed area. For Orion, the exposed area 


are the solar panels which are expected to be pointed toward the Sun. Therefore 


eg = 1s (50) 
Psp = 4.51 x 107° N/m? (51) 
A = 72.66 m? (52) 
m = 24080 kg (1650 slugs) (53) 
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This results in 

Fgpp = 1.105 x 107+ lbp (54) 
with the resulting acceleration 

asrp = 6.697 x 10-8 ft/s? (55) 


For a LINCOV integration step size of 60 seconds, the process noise is chosen to be (spreading it 


out in 3 axes) 


1 2 
Qsrp = (<zasne At = 8.970 x 107” ft?/s° (56) 


The Active Process Noise 


Finally, there is the case of the ammonia sublimator running occasionally. For EM-2 and fol- 
lowing, it is not expected that there will be any sublimator vents while in lunar orbit or close to 
the Moon. But the impulse prior to EI is expected to be 505.4 lbf-s for 20 minutes. The AV is 


expected to be 


FAt S 505.4 lbf — sec 


A i a 
ms m 673 slugs 


= 0.751 ft/s = 0.2289 m/s (57) 


This means that 


0.751 


2 
Q sublimecat onty — ( A ) /1200 sec = 1.5667 x 10~* ft?/sec®? = 1.4555 x 107° m?/sec3(58) 
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NAVIGATION AND DISPERSION RESULTS 
The Translunar Phase 


The navigation errors and trajectory dispersions for the EM-1 mission with optical navigation 
mapped to the B-plane are presented in Figures 3-4. For a sense of scale, these navigation errors 
and trajectory dispersions mapped to the B-plane with the (radius of the) Moon plotted for scale 


are presented in Figures 5-6. 
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Figure 3. The Translunar Navigation Errors mapped to the B-plane (with optical navigation) 


Finally, the performance at the flyby are presented in Figures 7 - 8. 


Type Lunar Flyby 30 Value 
Periapsis Altitude 0.675 n.m. 
Inclination 0.0238 deg 
Argument of Periapsis 0.0258 deg 
Cs 121890 ft? /s? 


Table 3. The Trans-Lunar Delivery 3c Statistics at OPF (with optical navigation) 
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Figure 4. The Translunar Trajectory Dispersion mapped to the B-plane (with optical navigation) 


Type Lunar Flyby 30 Value 

B-T 3.90 n.m. 

B-R 21.50 n.m. 
Linearized Time-of-Flight 0.814 sec 


Table 4. The Trans-Lunar Delivery 3c B-plane Statistics at OPF (with optical navigation) 


The TransEarth Phase 


The navigation errors mapped to Entry Interface flight path angle errors are presented in Figures 
9. The trajectory dispersions mapped to Entry Interface flight path angle trajectory dispersions 
are presented in Figures 10. Figures 11-13 contain the trajectory dispersion ellipses mapped to 
EI so that the correlations between the components (Downrange Position and Flight Path Angle, 
Velocity Magnitude and Flight Path Angle, and Crossrange Position and Crossrange Velocity) are 
observedThese plots also have the EI requirement displayed. One can see that that the requirement 
for Velocity Magnitude and Flight Path Angle is violated. The major contributor to this violation 


is the late ammonia sublimator vent that occurs 0.5 hour before EI. 
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Figure 5. The Translunar Navigation Errors mapped to the B-plane (with Moon for scale) 
(with optical navigation) 


CONCLUSIONS 

This paper has described the navigation and dispersion analysis performed to date on the Orion 
EM-2 mission, the first crewed mission. The analysis shows that the mission meets the Entry 
Interface Requirements, albeit barely. The models used are conservative and will be validated on 


EM-1. 
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Figure 6. The Translunar Trajectory Dispersion mapped to the B-plane (with Moon for 
scale)(with optical navigation) 
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Figure 7. The Translunar Trajectory Delivery Dispersion mapped to the B-plane (with optical 
navigation) 
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Figure 8. The Translunar Trajectory Delivery Dispersion mapped to the B-plane (with optical 
navigation) 
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Figure 9. The on-board Navigation Errors Mapped to Entry Interface Flight Path Angle Errors 
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Orion On-Board Trans-Earth Trajectory Dispersions Mapped to Entry Flight Path Angle 
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Figure 10. The Trans-Earth Trajectory Dispersions Mapped to Entry Interface Flight Path 
Angle Dispersions 
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Figure 11. The Trans-Earth Entry Interface Conditions Mapped Downrange vs Flight Path 
Angle Trajectory Dispersion Ellipses 
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Figure 12. The Trans-Earth Entry Interface Conditions Mapped Velocity Magnitude vs Flight 
Path Angle Trajectory Dispersion Ellipses 
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Figure 13. The Trans-Earth Entry Interface Conditions Mapped Crossrange Positon vs 
Crossrange Velocity Trajectory Dispersion Ellipses 
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Figure 14. The Trans-Earth Entry Interface Conditions Mapped Downrange vs Flight Path 
Angle Delivery Error Ellipses 
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Figure 15. The Trans-Earth Eniry Interface Conditions Mapped Velocity Magnitude vs Flight 
Path Angle Delivery Error Ellipses 
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Figure 16. The Trans-Earth Entry Interface Conditions Mapped Velocity Magnitude vs Down- 
range Delivery Error Ellipses 
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Figure 17. The Trans-Earth Entry Interface Conditions Mapped Crossrange Positon vs 
Crossrange Velocity Delivery Error Ellipses 
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